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We describe a simple, practical scheme for generating multi-qubit W states in resonator-based 
architectures, in which A'^ Josephson phase qubits are capacitively coupled to a common resonator 
bus. The entire control sequence consists of three pulses: a local Rabi pulse that excites a single 
qubit in the circuit; a coupling pulse that transfers the qubit excitation to the resonator bus; and 
the main, entangling operation that simultaneously couples the bus to all A'' qubits. If the qubit- 
resonator coupling strength g is much smaller than the qubit energy splitting _Eio, the system 
initially excited into the near-degenerate single-excitation subspace stays within that subspace, 
while smoothly evolving toward the fully uniform W state superposition. The duration of the 
final entangling operation is found to decrease with the total number of the qubits according to 
t = n/{2g\'N), in agreement with some of the previously proposed cavity QED W state generation 
schemes. 
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I. DESCRIPTION OF THE W PROTOCOL 

Our control sequence consists of the follovifing three 
steps: 



T 



First, the initial local Rabi pulse is applied to one 
of the qubits in the circuit, bringing the qubit from 
its ground state |0) to the excited state |1), 



|00...000^) -^ |00...010^). 



(1) 



2. The corresponding qubit-resonator coupling is then 
turned on, which transfers the excitation state from 
the qubit to the bus, 



|00...010^) -^ I00...0010. 



(2) 



The second entangling pulse is applied, which cou- 
ples the bus to all the qubits in the system. If the 
coupling g is much smaller than the qubit energy 
splitting EiQ, the system initially prepared in the 
single-excitation subspace stays within that sub- 
space while smoothly evolving toward the multi- 
qubit, fully uniform W-state superposition as fol- 
lows. 



00... 001,.) -^ 

" 00... 01) + 100...10) + ••• + 110...00) 



TV 



«) 10^ 



(3) 



This last step is similar to the W state generation 
scheme proposed for cavity QED in Ref. p,]. In 
the terminology of reference [2| , our resonator bus 
plays the role of the entanglement mediator. 
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II. MATHEMATICAL PRELIMINARY 

It is well-known how to perform the first two opera- 
tions of the W sequence described above [3]. We can 
therefore assume that the circuit was initially prepared 
in the state |00 . . .001^), with only the bus excited. By 
simultaneously turning on the N couplings, the W state 
of the A^-qubit network can then be generated using a 
single entangling operation (cf. [1]). 

In order to see how this works, we consider a formal 
problem of an "effective" Hamiltonian, 
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(4) 



which operates within a certain {N + l)-diniensional 
Hilbert space 'H'-^^^\ whose physical significance will 
be clarified below. The spectrum of if (^+i) is found to 
be 



ij(^+i) = ^VTv,o,...,o. 



(5) 



The corresponding eigenvector matrix S^ ' , which di- 
agonalizes i/'^+D via HJ,ZV'> = ^f^+Dti/C^+D^^^+D, 
is given by 
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(6) 



where we left the columns of 5'^^"'"^^ unnormalized for 
notational simplicity. Direct exponentiation then shows 



that the TV-dimensional uniform superposition state in 
this "effective" {N + l)-dimensional system can be gen- 
erated via 
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-iff("+i)t<"' 
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(7) 



where 



t 



(N) ^ 



2gy/N 



(8) 



III. Wjv STATE GENERATION 

Our W state generation scheme is based on the 
idea that the effective Hamiltonian iy(^+i) considered 
above should be viewed as operating within the single- 
excitation subspace of a network consisting of N qubits 
coupled to a common resonator bus. The corresponding 
mapping (extended by linearity) between the "effective" 
Hilbert space 7^(^+1) and the (A^-l-l)-dimensional single- 
excitation subspace of the system may be chosen to be 
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I00...001, 
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|00...010^ 



(9) 



The uniform A^-dimensional superposition state gener- 
ated in accordance with Eq. (|22|) then corresponds to the 
Wat state of the N qubits attached to the common bus. 

In order for this approach to work, the single-excitation 
subspace has to be well isolated from the rest of the sys- 
tem's Hilbert space. This near-degeneracy condition is 
typically well satisfied in various superconducting qubit 
architectures whose couplings, g ~ 100 MHz, are much 
smaller than the qubit and resonator level splittings of 
Eio ^ 10 GHz. 

Let us check that the Hamiltonian iy(^+i) arises nat- 
urally within the single-excitation subspace of a capac- 
itively coupled network consisting of superconducting 
phase qubits and a resonator bus. When projected into 
the computational subspace spanned by the eigenfunc- 
tions |0), |1), |2) of the individual Josephson phase qubits 
(as well as the resonator) , the Hamiltonian of such a net- 
work is given by 



N 



H = Y,H.^ + Hr 



N 



dirPiPr 



(10) 



where the index i numbers the qubits and r labels the 
bus, with 



Hi = 



Hj = 



Hr — 



'^10 



^10-Ai, 



-E 




J =2, 



A, 



,iV, 



(11) 



The generalized momenta pi and Pr are given by 



/O -1 -&, 

Pi = A2 + bXr-, + cXj = i 1 ~Ci 

\bi Ci 

/O -1 

Pr = X2+ \/2A7 = i 1 -V2 
VO %/2 



(12) 



where \k, k ~ 1,2,. ..,8, are the standard Gell-Mann 
generators of the Lie algebra su(3). In the above, gir are 
the qubit-bus coupling constants, Eiq is the energy split- 
ting of the first (reference) qubit, e^, e^, j = 2, 3, . . . , A^, 
are the energy shifts relative to the single-excitation en- 
ergy of the first qubit, A^, i = 1,2,..., TV, are the 
qubit anharmonicities, bi and Ci are the off-diagonal ma- 
trix elements of the ith qubit momentum, and £",- is 
the resonator energy splitting. The (TV -f 1) x (N + 1) 
block of the Hamiltonian H acting within the {N + 
l)-dimensional single-excitation subspace spanned by 
|00... 001^), |00... 0100,..., |10...000r), is then given 
by the real symmetric matrix. 
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/ Cr gNr gN-lr gN~2r 

gNr CAT 

gN-lr EAT-I 

gN-2r eAr-2 



53r 52r fflrX 









i=l 



52r ... £2 

\ gir ...000/ 

(13) 
This immediately shows that the Wat state can be gen- 
erated if we place all system elements on resonance with 
each other by choosing 



£2 = £3 = £4 = • • • = ew = Cr = 0, 

gir = 52r = 53r = • • • = gNr = 5- (14) 



IV. NUMERICAL SIMULATION RESULTS 

We tested this scheme on an iV = 4 qubit network with 
g = 100 MHz, £'10 = Er = 10 GHz, A^ = 250 MHz, and 



-1 -0.08^ 
Pi=i \ 1 -1.43 
,0.08 1.43 



(15) 



Assuming the system starts in the excited state 
|00. . .001,-), the simulated final state of the system is 
found to be 



|WAr)sim = 



/-O.OOOSA 

' 0.4999 ' 

0.4999 

0.4999 

V 0.4999 / 



(16) 



with the corresponding entangling time being t = 1.2500 
ns. Ignoring the decoherence effects, the intrinsic fidelity 
[J| of the found state |WOv)sim relative to the ideal Wat 
state, is 



.F^ = |(WAr|WAr),i„,|2 = 0.9994. 



Wat+i state generation 



In a similar manner, the Wat+i state can also be gener- 
ated, in which the resonator is maximally entangled with 
the qubits. This corresponds to the sequence of opera- 
tions 

I00...000,.) ^ |00...010^) -^ 

|00 . . . 001^) -h |00 . . . 010,.) + • • ■ + |10 . . . 000^) 



^/N+l 



(18) 



In this scenario, we take full advantage of qubit tunability 
to construct the "effective" single-excitation Hamiltonian 
of the form 
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Its spectrum and the diagonalizing transformation ma- 
trix S'(^+i) (here shown unnormalized) are 
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E("+') = l=FVAf + l,0,...,0, 
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(21) 
respectively. The corresponding W state is generated via 
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(17) where 



JN+l) _ 



^(W+1) = 



(22) 
(23) 



2y/N + 1 2g^/NTT 

provided we detune the qubits from the resonator by 2g 



£2 = £3 = £4 = 

fflr = 32r = .93r 



En = 0, €,. = 2.9, 

■ = gNr = .9- 



(24) 



The duration of the entangling pulse is now t = 1.1180 
ns, with the simulated single-excitation final state of the 
network being 
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with fidelity 
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